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ABSTRACT 

A problem still unsolved in cosmology is the identification of the sources of radiation 
able to reionize H I in the intergalactic medium (IGM) by z ~ 6. Theoretical works 
and observations seem to indicate that the fraction, (f esc ), of Hi ionizing radiation 
emitted from galaxies that escapes into the IGM is small in the local universe ((fesc) 
£ 10%). At high redshift galaxies are more compact and probably gas rich implying 
smaller values of (fesc) from their disks or spheroids. But if the sites of star formation 
are displaced from the disk or spheroid and the star formation efficiency of the proto- 
clusters is high, then (f esc ) should be about one. This star formation scenario is 
consistent with several models for globular clusters formation. Using simple arguments 
based on the observed number of globular cluster systems in the local universe and 
assuming that the oldest globular clusters formed before reionization and had (fesc) 
~ 1, I show that they produced enough ionizing photons to reionize the IGM at z w 6. 
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1 INTRODUCTION AND RATIONALE 



Observation of Lya ab sorption systems toward newly found 



high- redshift quasars (Becker et al. 2001; Djorgovski et al 



|200l| ) indicate that the redshift of reionization of the inter- 



galac t ic medium (IGM) shou ld be close to z = 6 (Gnedin 
2002; Songaila & Cowie 2002). Perhaps the recent identifi- 



large: (/esc) ~ 10% assuming a Salpeter initial mass func- 
tion (IMF) and the standard ACDM cosmological model. 
Benson et al. (2002) finds that (fesc) should be about 15% 
for reionization at z = 6, but a smaller value (f esc ) ~ 10% is 
consistent wit h the observed ionizing background at z ~ 3. 



Gnedin (2002) finds that assuming a primordial power spec- 



cation of a lensed galaxy at z = 6.56 points to a somewhat 
earlier redshift of reionization ( Hu et al. 2002 ). Although 
quasar s play a dominant role in photoionizing the IGM at 
z ~ 3 (Meiksin & Madau 1993), their dwindling numbers at 
z > 4 suggest the need for another ionization source. Unless 
a hidden population of quasars is found, radiation emitted 
by high-redshift massive stars seems necessary to reionize 
the universe. A key ingredient in determining the effective- 
ness by which galaxies photoionize the surrounding IGM is 
the parameter (fesc), defined here as the mean fraction of 
Lyc photons escaping from galaxy halos into the IGM. To 
be an important source of ionizing photons and rival with 
quasars, a substantial fraction (~ 10%) of them mus t escape 
the gas layers of the galaxies (Madau & Shull 1996). 

Cosmological simulations and semi-analytical models 
of IGM reionization by stellar sources find that the ioniz- 
ing background rises steeply at the redshift of reionization. 

TJ nf ort nnately a direct pnmparisnn het.wppn models is dif- 



trum index n = 0.93, the preferred value from CMB and 
LSS data, reionization at z £ 6 requires a large (fesc)', but 
this assumption produces an ionizing background at z £ 4, 
that is too large. The common assumption of a universal 
star formation efficiency (SFE) (e.g., the coefficient in front 
of the Schmidt-Law in some models or in others the fraction 
/» of baryons converted into stars) is consistent with the ob- 
served values of the star formation rate (SFR) at < z < 5 
and total star fraction f2„ at z — 0. However, the assump- 
tion of a constant (fesc) does not seem to be consistent with 
observations. An escape fraction (fesc) ~ 1 is required for 
reionization at z ~ 6 but the io nizing background a t z ~ 3 
is consistent with (fesc) ~ 10% (Bianchi et al. 2001). Small 



values of (fesc) at z £ 3 are also supported by direct obser- 
vations o f the Lyc emission from Lyman-break and starburst 
galaxies. Giallongo et al. (2002) find an upper limit (fesc) 
< 16% at z rZ 3 (but see steidel et al. 2001) and obser- 
vations of low-redshift starbursts arc consistent wi th (fesc) 



upper limits ranging from a few pe rcent up to 10% ( Hurwitz 



ficult b lecause of different recipes used for star formation. 
clumping of the IGM or the definition of (fesc). But a result 
common to all the models is that, in order to reionize the 
IGM by z = 6 — 7, the escape fraction must be relatively 



et al. 1997; Deharveng et al. 2001) 



Calculations of (fesc) from first principles are difficult. 
The main complications arise in simulating a realistic in- 
terstellar medium (ISM) that includes small-scale physics 
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and feedback processes. Moreover the mean (fesc) results 
from the contribution of a variety of galaxies whose ISM 
properties ar e largely unknow n at high redshift. Theoret- 



ical models (Dove et al. 2000 CJiardi et al. 2002) for the 



radiative transfer of ionizing radiation through the disk 
layer of spiral galaxies similar to the Milky Way find (/esc) 
~ 6 — 10%. At high redshift the mean value of (fesc) is ex- 
pecte d to decrease almost e x ponentially with incr easing red- 
shift ( |Ricotti fc Shull 2000| ; [Wood fc Loeb 2000| ); at z > 6, 
(fesc) ~ 0.1 — 1% even assuming star formation rates typ- 
ical of starburst galaxies (e.g., SFR~ 10 times that of the 



Milky Way). Using Monte-Carlo simulations, Ricotti & Shull 
(200o|)~[iave studied how (fesc) depends on galactic param- 
eters. Assuming gas density profile in hydrostatic equilib- 
rium in the dark matter (DM) potential, star density pro- 
portional to the gas density and a power law for the lu- 
minosity function of the OB associations, they found that 
(fesc) cx (ef g M D M)~ 1/3 exp[-(z vir + l)e" 1/3 ]. Here e is pro- 
portional to the SFE, f g is the fraction of collapsed gas, z V ir 
is the virialization redshift and Mdm is the DM halo mass. 
The majority of photons that escape the halo come from 
the most luminous OB associations located in the outer- 



most parts of the galaxy. Indeed Ricotti & Shull (2000) have 
shown that changing the luminosity function of the OB as- 
sociation and the density distribution of the stars has major 
effects on {fesc) (see their Figs. 8 and 9). In the aforemen- 
tioned models, (fesc) should be regarded as an upper limit 
since dust extinction and absorption of ionizing radiation 
from the molecular cloud in which OB associations are born 
are neglected. 

The theoretical suggestion of a decreasing (fesc) with in- 
creasing redshift is in contrast with models for reionization 
that require (fesc) ~ 1 at z — 6. A different star formation 
mode, with very luminous OB associations forming in the 
outer parts of galaxy halos, could explain the large (fesc) re- 
quired for reionization. Globular clusters (GCs) are possible 
observable relics of such a star formation mode. Their red- 
s hift of formation is compatible with redshift of reionization 
( Gnedin et al. 2002 ) . Because of their large star density they 
survived tidal destruction and represent the most luminous 
tail of the luminosity distribution of primordial OB associa- 
tions. In § 2.1 I explain that several models for the formation 
of proto-GCs imply an (fesc) ~ 1. I will also show that the 
total amount of stars in GCs observed today is sufficient to 
reionize the universe at z ~ 6 if their (fesc) ~ 1. This con- 
clusion is reinforced if the GCs we observe today are only a 
fraction, 1/fdi, of primordial GCs as a consequence of mass 
segregation and tidal stripping. 

The paper is organized as follows. In § |^ I briefly review 
recent progress in our understanding of GC properties and 
formation theories; in § ^ I discuss the model assumptions 
in light of GC observations and present the results. In § ^ I 
present my conclusions. 



2.1 Observations 

Most galaxies have a bimodal GCs distribution indicat- 
ing that luminous galaxies experience at least two major 
episodes of GCs formation. The bulk of the globulars in the 
main body of the Galactic halo appear to have formed during 
a short-lived burst (~ 0.5 — 2 Gyr) that took place about 
13 Gyrs ago. This was followed by a second burst associ- 
ated with the formation of the galactic bulges. Clusters may 
have been formed in d warf spheroidal galax ies and acreted 
by the Galactic halo ( van den Bergh 1999| ). Massive clus- 
ter formation occurred in galaxies as small as the Fornax 
dwarf spheroidal but not in massi ve ones such as the small 
magellanic cloud ( Zepf et al. 1999] ). 



2.1.1 Absolute and relative ages 

The method for determinating the absolute age of GCs is 
based on fitting the observed color-magnitude diagram with 
theoretical evolutionary tracks. The systematics in the evo- 
lutionary model and the determination of the cluster dis- 
tance are the major sources of errors. Recent determinations 
of the absolute age of old GCs find t gc = 12.5 ± 1.2 Gyr 
( |Chaboyer et al. 1998 ) , consistent with radioactive dating 
of a very m etal-poor star in the halo of our galaxy (Cayrel 
et al. 2001). Relative ages of Galactic GCs can be deter- 



mined with greater accuracy, since many systematic errors 
can be eliminated. In our Galaxy At gc = 0.5 Gyr, but differ- 
ences in age betw een GC systems in d ifferent galaxies could 



be At g c ~ 2 Gyr ( |gtetson et al. 199fi| ). 



2.1.2 Specific frequency 

The specific frequency is defined as the number, N, of 
GCs per My = —15 of parent galaxy light, Sn = N X 
10 o.4(jvv+i5) ( Harris & van den Bergh 1981). The most 

striking characteristic is that Sat (Ellipticals) > Sat (Spirals). 
Sn = 0.5 in Sc/Ir galaxies (Harris, 1991), Sn = 1 in spi- 
rals of types Sa/Sb , Sn = 2.5 in field ellipticals (Kundu 



fc Whitmore 2001e] , bj). Converting to luminosity (LvJLq = 
10 -u.4im„-4. MJ j ^ e have N = (L/L Q )S N /8.55 x 10 7 . I can 
therefore calculate the efficiency of GC formation defined as, 



Mgc 
M 



fdiNm g 
M 



SNfdi 

(M/L) v 



x 0.00585, 



(1) 



where M gc is the total mass of the GC system, m ffc = 5 x 
10 Mq is the mean mass of GCs today, M is the stellar 
mass and (M/L)v is the mass to light ratio of the galaxy. 
In the next paragraph we show that, because of dynamical 
evolution, m gc and iV are expected to be larger at the time of 
GC formation than today. Therefore, the parameter fdi ^ 1 
is introduced to account for dynamical disruption of GCs 
during their lifetime. 



2 SHORT REVIEW ON GC SYSTEMS 

In this section I review some observational and theoretical 
results on GCs useful to the aim of this paper. I also try 
to justify my assumption (fesc) ~ 1 for GCs on the base of 
theoretical models of proto-GC formation. 



2.1.3 IMF, Metallicity and Dynamical Evolution 

The IMF of GCs is not known. The present mass function is 
known only between 0.2 and 0.8 Mq since high-mass stars 
are lost because of two-body relaxation and stellar evolu- 
tion processes. Theoretical models show that the shape is 
consistent with a Salpeter-like IMF. The mean metallicity 
of old GCs is Z ~ 0.03 Zq. One of the most remarkable 
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properties of GCs is the uniformity of their internal metal- 
licity A[Fe/H] £0.1. This implies that the bulk of the stars 
that constitute a GC formed in a single monolithic burst of 
star formation. A typical GC emits S ~ 3 x 10 53 s _1 ioniz- 
ing photons in a burst lasting 4 Myrs: about 300 times the 
ionizing luminosity of largest OB associations in our Galaxy. 

During their lifetime GCs lose a large part of their ini- 
tial mass or are completely destroyed by internal and ex- 
ternal processes. Prominent internal processes are mass loss 
by stellar evolution (about half of the initial mass is lost) 
and two-body relaxation, the effects of which are mass seg- 
regation (change in the mass function) and core collapse 
(expansion and evaporation). External processes can be di- 
vided roughly into two classes: gravitational shocks from GC 
motion through the disk or bulge of the galaxy, and tidal 
forces, which cause mass loss due to tidal truncation. Nu- 
merical simulations show that about 50%-90% of the mass 
of GCs is lost due to external processes, depending on the 
host galaxy env ironment, initial concentrat i on and IMF of 



Table 1. Star census at z=0. 



the proto-GCs (Chernoff fc Weinberg 199C; Gmedin fc Os- 



triker 1997). Many of the low-metallicity halo field stars in 
the Milky- Way could be debris of disrupted GCs. The mass 
in stars in the halo is about 100 times the m ass in GCs. 
Therefore the parameter fdi, defined in § 2.1.2, could be as 



large as fdi = 100. Overall fdi is not well constrained since it 
depends on unknown properties of the proto-GCs. Accord- 
ing to results of N-body simulations fdi should be in the 
fdi ~ 2 - 10 range. 



2.2 Why is (f e . 



1 plausible for GCs? 



I discuss separately two issues: (i) the (fesc) from the gas 
cloud in which the GC forms, and (ii) the (/esc) through any 
surrounding gas in the galaxy. 

(i) The evidence for (fesc) ~ 1 comes from the observed 
properties of present-day GCs. The fact that they are com- 
pact self-gravitating systems with low and uniform metal- 
licity points to a high efficiency of conversion of gas into 
stars. A longer timescale of star formation would have en- 
riched the gas of metals and the mechanical feedback from 
SN explosions would have stopped further star formation. 
If /* w 10% of the gas is converted into stars in a single 
burst (with duration < 4 Myr) at the center of a spheroidal 
galaxy, follo wing the simple calculations shown in Ricotti & 
Shull (2000| ) [see their eq. (18)] at z = 6 we have, 



fesc = 1 ~ 0.06 



7 



50%. 



(2) 



(ii) The justification for (fesc) ~ 1 is model dependent but 
in general there are two main arguments: a) the high effi- 
ciency of star formation /„, and b) the sites of proto-GC 
formation in the outermost parts of the galaxy halo. 
In the "cosmo logi cal objects model" (30 < Zf < 7) of 



Pcc- 



bles fe Dicke (1968) GCs form with efficiency /» w 100%, 



implying (fesc) = 1 (note that such a high /* is not found in 



numerical simulations of first object formation ( Ricotti et al 



2002a 


b|)). In "hierarchical formation models" (10 < Zf < 3) 


(ILarsc 


m 1993; Harris & Pudritz 1994 




McLaughlin & Pudritz 


199f ; Gnedin et al. 2002 ) GCs form in the disk or spheroid 



of galaxies with gas mass M g ~ 10 ' - 10 a M . Compact 
GCs survive the accretion by larger galaxies while the rest 



Type 


u* (%) 


Sn 


(M/L) v 


£ S c (%) 


Sph 


6.5+ 4 3 


2.4 ±0.4 


5.4 ±0.3 


0.26 ± 0.06 


Disk 


n+1.5 
z -0.5 


1±0.1 


1.82 ±0.4 


0.32 ± 0.1 


Irr 


15+ 015 
U - L0 _0.05 


0.5 


1.33 ±0.25 


0.22 ± 0.04 


Total 


Q + b.b 

a -3.5 






0.3 ±0.07 



of the galaxy is tidally stripped. Assuming that 1 — 10 GCs 
form in a galaxy with M g ~ 10 7 - 10 8 M implies /* ~ 10% 
and therefore (fesc) ^ 50%. (fesc) is larger than in eq. ^ 
if proto-GCs are located off-center (e.g., if they form from 
cloud-cloud collisions during the galaxy assembly) or if part 
of the gas in the halo is collisionally ionized as a consequence 
of the virialization process. 

In models such as the "s upershell fragmentation" (zf < 10) 
of Taniguchi et al. (1999 ) or the "thermal instability" (zj < 
7) of [Fall fc Rees (1985fl ! (fesc) ~ 1 since proto-GCs form in 



the outermost part of an already collisionally ionized halo. 

In summary, since (fesc) depends strongly on the luminosity 
of the OB associations and on their location, proto-GCs, 
being several hundred times more luminous than Galactic 
OB associations, should have a comparably larger (fesc). 



3 METHOD AND RESULTS 

In this section I estimate the number of ionizing photons 
emitted per baryon per Hubble time, JV p h, by GC formation. 
In § [0| I derive Af p h assuming that all GCs observed at z = 
form ed i n a time period At gc with constant formation rate. 
In § |3.2| I use the Press-Schechter formalism to model more 
realistically the formation rate of old GCs. 



3.1 The simplest estimate 

I start by estimating the fraction, uj gc , of cosmic baryons 
converted into GC stars. By definition ui gc = uj,e gc , where 
w* is the fraction of baryons in stars at z — and e gc is 
the efficiency of GC formation defined inj 2.1.2. In all the 
calculations I assume Qb = 0.04. I n Tabic pL I summariz e the 
star census at z — according to Fukugita et al. (1998) and 
I derive e a c using eq. (^), ass uming fdi = 1. Using similar 
arguments McLaughlin (1999) finds a universal efficiency of 
globular cluster formation e gc = (0.26±0.05)%, in agreement 
with the simpler estimate presented here. It follows that 
+l 3 7 x 10- 4 ^ 



gc = /<k(2.7I£? x KT 4 ) at z = 0. 

The total number of ionizing photons per unit time 
emitted by GCs is rjuj g e/ At gc , where r\ is the number of ion- 
izing photons emitted per baryon converted into stars, and 
(jj gc ~ 2.1w 9C takes into account the mass loss due to stellar 
winds and SN explosions adopting an instantaneous-burst 
star formation law. GCs did not recycle this lost mass since 
they formed in a single burst of star formation, n depends 
on the IMF and on the metallicity of the star. I calculate n 
u sing a Salpeter IMF and star metallicity Z — 0.03 Z q (see 
§ 2.1.3) with Starburst99 code ( Leitherer et al. 1999] ), and 
find rj = 8967. The number of ionizing photons per baryon 
emitted in a Hubble time at z — 6 is, 



tjajz = 6) 

Atac 



(5.1 



+4.3 n 
3.2/ 



fdi 



At g e (Gyr) ' 



(3) 
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where I have assumed (fesc) = 1 and Hubble time at 2 = 6 
t H = 1±0.1 Gyr. I expect 1 < f di < 100 and 0.5 < At gc < 2 
Gyr. A conservative estimate of fdi ^. 2 and At gc £ 2 Gyr 
(i.e., 10 < zf < 3) implies f dt /At gc > 1. The IGM is reion- 
ized when A/" p h = C, where C = {n 2 HII ) / (uhii) 2 is the ion- 
ized IGM clumping factor. According to the adopted defini- 
tion of (fesc), C = 1 for a homogeneous IGM, or 1 ,$ C ^ 10 
taking into account IGM density fluctuations producin g the 
Lya forest ( |Miralda-Escude et al. 200C| ; |Gnedin 200o| ). The 
estimate from eq. (^) is rather rough because I have implic- 
itly assumed that the SFR is constant during the period of 
GC formation At gc . A more realistic SFR as a function of 
redshift requires assuming a specific model for the formation 
of GCs. I try to address this question in the next section. 



3.2 A bit of modeling 

I assume that the formation rate of stars or GCs in galaxies is 
proportional to the merger rate of galaxy halos (each galaxy 
undergoes a major star burst episode when it virializes). 
Using the Press-Schechter formalism I calculate, 



10.23 



t (Gyr) 
12.54 13 



13.16 



du s gc (z) 
dt 

dudjz) 
dt 



A 



B 



M2 dQ{M DM ,z) 
dt 

dSl(M DM ,z) 
dt 



Mi 



din M dm, 



din Mdh 



(4) 



(5) 



where Q(Mdm , z)dln(MoM) is the mass fraction in viri- 
alized DM halos of mass Mdm at redshift z. I determine 
the constants A and by integrating eqs. (Q)-(^|) with re- 
spect to time, and assuming uj gc — 0.1% (i.e., fdi = 2) and 
loI — 1.4w* = 13% at 2 = (the factor 1.4 takes into account 
the mass loss due to stellar winds and SN explosions adopt- 
ing a continuous star formation law). I assume that GCs 
form in halos with masses M\ < Mdm < M2- The choice of 
Mi and M2 determine the mean redshift, 2/, and time pe- 
riod, At 9 c, for the formation of old GCs. In order to be con- 
sistent with observations I consider three cases: case (i) halos 
with virial temperature 2 x 10 4 < T V i r < 5 x 10 4 K; case (ii) 
5x 10 4 < T vir < 10 5 K; and case (iii) 10 5 < T vir < 5x 10 5 K. 
In case (i) , (ii) and (iii) At gc — 2.2, 3.7 and 5.2 Gyr, respec- 
tively, and the GC formation rate has a peak at z = 7.5, 6 
and 4.6, respectively. Disk and spheroid stars form in halos 
with Mdm > M m . At 2 > 10 I assume that the first ob- 
jects form in halos with M m corresponding to a halo virial 
temperature T V i r = 5 x 10 3 K. A t 2 < 10 only object s with 
T vir > 2 x 10 4 K can form (see [Ricotti et al. 2002b| ). The 
comoving star formation rate, given by p„ = ~puj{ , where 
7j = 5.51 x 10 9 M Q Mpc" 3 is the mean baryon density at 
2 = , is shown in Fig, hi Th e points show the observed SFR 
from Lanzetta et al. (2002). In Fig. I show M v h for GCs 



(thick lines) and for galaxies (thin lines) defined as, 



, dw s gc 

nfdi—^-t H (z), 



dt 



X 



J 



P h 



V(fesc)^-t H {z). 



(6) 

(7) 



1 By definition J °° U(M DM , z)<1\xi(M d m) = 1. I find the follow- 
ing values of the constants A = 1.3%, 1.6%, 0.6% for cases (i), (ii) 
and (iii) respectively (see text) and B = 12%. 




Figure 1. The solid line in the top panel shows the comoving 
SFR as a function of time in our model. The solid, dashed and 
short-dashed lines show the SFR of GCs for cases (i), (ii) and 
(iii), respectively (assuming fdi = 2). The bottom panel shows 
as a function of time. The segment with arrows is a visual aid 



to compare oj gc (assuming fdi = 20) to around z ~ 6 



The thick solid, dashed and short-dashed lines show for 
case (i), (ii) and (iii), respectively. For comparison, I show 
(thin solid line) J\f* h assuming (fesc) = 0.1 x exp[— 2/2], 
chosen to fit the observed values ( squares) of Af£h at z = 
2,3,4 ( |Miralda-Escude et al. 200C ). The thin dashed line 
shows M* h assuming constant (/ esc ) = 5%. 



4 CONCLUSIONS 

Observed Lyman break galaxies at 2 ~ 3 are probably the 
most luminous starburst galaxies of a population that pro- 
duced the bulk of the stars in our universe. Their forma- 
tion epoch corresponds to the assembly of the bulges of spi- 
rals and ellipticals. Nevertheless the observed upper limit 
on (fesc) from Lyman break galaxies is (fesc) ~ 10%, insuf- 
ficient to reioni ze the IGM according to numerical simula- 
tions. Recently Ferguson et al. (2002), using different argu- 
ments, have claimed that the radiation emitted from Lyman 
break galaxies is insufficient to reionize the IGM assuming 
a continuous star formation mode. 

I propose that GCs could produce enough ionizing pho- 
tons to reionize the IGM. Assuming fdi = 2 (i.e., during 
their evolution GCs have lost half of their original mass), I 
find w£ c » 0.1%, small compared to the total u>{ ~ 10% at 
2 = 0. But GCs are around 12-13 Gyr old and, if they formed 
between 7 < 2 < 5 (in about 0.5 Gyr), the expected total oj( 
formed during this time period is about ul ~ 1%, only 10 
times larger than u> gc . Assuming fdi = 20, expected from the 
results of N-body simulations, I find cu gc ~ 1%, suggesting 
that GC formation is an important mode of star formation 
at high-redshift. The special star formation mode required 
to explain the formation of GCs suggests an (fesc) ~ 1 
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2 4 6 8 10 12 14 



z 

Figure 2. Emissivity (photons per baryon per Hubble time) as 
a function of redshift. The thick solid, dashed and short-dashed 
lines show the contribution of GCs for cases (i), (ii) and (iii), 
respectively (assuming = 2). The thin lines show the contri- 
bution of galaxies assuming a realistic {fesc) = 0.1 X exp[— z/2] 
(solid) and a constant (fesc) = 5% (dashed). 

from these objects. This is because {f e3C ) is dominated by 
the most luminous OB associations and GCs are extremely 
luminous, emitting S ~ 3 x 10 53 s _1 ionizing photons in 
bursts lasting 4 Myrs. Moreover, according to many mod- 
els, GCs form in the hot, collisionally-ionized galaxy halo, 
from which all the ionizing radiation emitted can escape into 
the IGM. Therefore it is not too surprising, if GCs started 
forming before z — 6, that their contribution to reioniza- 
tion is large. I find that the number of ionizing photons 
per baryon emitted in a Hubble time at z — 6 by GCs is 

= (S.li^'D/di/Atgc > 1, therefore sufficient to reionize 
the IGM even if we assume fdi = 1. Here, At gc ~ 0.5 — 2 is 
the period of formation of the bulk of old GCs in Gyrs. Using 
simple calculations based on Press-Schechter formalism (see 
Fig. ^) I find that, if normal star formation in galaxies have 
(fesc) ~ 5%, GC contribution to reionization should be im- 
portant. If GCs formed by thermal instability in the halo of 
T V i r ~ 10 s K galaxies (case (iii)), the ionizing sources are 
located in rare peaks of the initial density field. Therefore, 
the mean size of intergalactic H II regions before overlap is 
large and reionization is inhomogeneous on large scales. 

In this letter I have considered the possibility that an in- 
creasing (/esc) at z ~ 6 due to GCs formation could explain 
IGM reionization and still be consistent with the observed 
values of the ionizing background at z < 3. Alternatively an 
increasing production of ionizing photons per baryon con- 
verted into stars, due to a varying IMF, would have similar 
effects on the IGM. Chemical evolution studies should be 
able to distinguish between these two scenarios. 
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